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Abstract 

Purpose: To investigate the Lie groups of transformations (symmetries) approximately admitted by a 
finite difference discretization of a nonlinear advection equation on an unstructured mesh. 
Method: A recently developed technique of "approximate group analysis" is applied to a differential 
approximation (otherwise referred to as an equivalent differential equation) corresponding to the above 
finite difference approximation. 

in . 

Results and Conclusions: We determine which groups from the infinite variety of groups admitted by 
a nonlinear advection equation "survive" the discretization. The situations arising for different choices 
of an arbitrary function (local speed of propagation) are also studied. 



o 
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Introduction 



Invariance of differential equations with respect to one-parameter groups of transformations in the space 
of independent and dependent variables carries important information about the fundamental properties 
of the physical system that these equations describe, such as conservation laws. Knowledge of the groups 
admitted by the system of differential equations may also allow one to reduce the order of the system, to 
find important particular solutions, or to produce new solutions from a single solution which is known JjJ. 

Preserving the group properties of differential equations in their discretizations seems to be a natural 
requirement which would ensure that the fundamental properties of the approximated differential equations 
are preserved in the discrete models B H pi [7J . This idea is often followed in the design of approximations 
used in CFD by making sure that the scheme preserves steady-state or lower-dimensional solutions exactly. 
Further illustrating this statement, the author of || refers to the result from M where it was demonstrated 
that the undesirable process of self-oscillation observed in the solutions obtained using Harlow's scheme M 
is caused by the non-invariance of the scheme with respect to the Galileo transformation. Unfortunately, 
in general, replacing the system of differential equations with their finite difference, finite volume or finite 
element discretizations introduces the non-invariance into the discretization. This non-invariance stems 
from the fact that performing the discretization requires utilizing a discrete set of nodes in the space of 
independent variables, i.e. a computational grid. Formally, this implies that the group of transformations 
acting on the discrete variables needs to leave invariant not only the approximating discrete system, but also 
a set of algebraic equations defining the grid |7j . Although some of the important groups of transformations 
arising in applications do satisfy this property, generally speaking, the discrete system loses many of the 
groups of the original differential equations. 

The natural way to avoid this restrictiveness refered to above is to weaken the requirement of invariancy 
of the discrete system. One of the possibilities is to consider the differential approximation Q or equivalent 
hyperbolic differential equation (EHDE) || corresponding to the discrete system. By formally applying a 
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Taylor series expansion to the approximating discrete system, one can obtain a modified differential equation 
which can be treated as a perturbation of the original differential equation due to discretization. 

The purpose of this work is to extend the above approach using a recently developed concept of approx- 
imate invariance JTo| , combining the methods of group analysis and perturbation theory. 

2 Method 

2.1 Nonlinear advection equation and its symmetries 

We consider a nonlinear advection equation 

L(x,t,u,ut,u x ) = u t + <p(u)u x = 0, (1) 

where (p(u) is a given arbitrary smooth function, characterizing the local speed of propagation of small 
disturbances. 

Let z = (t,x,u) be a point in the space of independent variables and unknown functions. Consider the 
following one-parameter family of transformations acting on z: 

z' = f(z,a), (2) 

where f = ifiifiifs) and a is a scalar parameter of this family of transformations. We say that a trans- 
formation (y|) is a one-parameter group (o.p.g.) (or a symmetry) with respect to the parameter a if (i) 
f (z, a) = z for all z, if and only if a = 0, and (ii) f (f (z, a), b) = f (z, a + b) for all a, b and z. It is convenient 
to associate each o.p.g. with its infinitesimal operator 

X = £ T d z = e (t) (x, t, u)d t + £ (x) (x, t, u)d x + V (x 1 t,u)d u , 

where £ = (^ t \ rj) T = -§-f\ a —Q- Lie theorem ||] establishes a one-to-one correspondence between an 
o.p.g. and its infinitesimal operator. 

Once the o.p.g. of transformations (j^) in the space of independent variables and unknown functions 
is specified, one can easily derive the corresponding transformations of the partial derivatives of unknown 
functions with respect to the independent variables. For equation (^), this procedure applied to partial 
derivatives of first order, results in an o.p.g. of transformations acting in the space of variables (x, t, u, Ut, u x ), 
which is called the first prolongation of the o.p.g (Q). One of the conveniences of using the infinitesimal 
operator to describe the o.p.g. is that there exist a compact prolongation formulae [jj], Theorem 2.36], which 
gives the components of the first prolongation X^ of the infinitesimal operator in terms of the components 
of the operator X . For the case under consideration, these formulae give: 

jfW =x + c (t) d Ut +c (x) d Ux , 

where 

C (t) = D t ( V - fWttt - + ^u tt + ^W, C (x) = DAV - ^ {t) ut - & ] u x ) + £Wufa + ^u xx , 

and D t = d t + u t d u + u tt d Ut + u tx d Ux + D x = d x + u x d u + u tx d Ut + u xx d Ux + are the first total 
derivatives. 

The main problem of group analysis of differential equations is to find all o.p.g. with respect to which 
the given differential equation is invariant. The solution to this problem can generally be found using 
infinitesimal criterion of invariance p], Theorem 2.31] which, for the present case, reduces to the following 
statement. Equation ([!]) is invariant with respect to an o.p.g. (^) if and only if 

xW(u t + ip(u)u x )\ Ut ^ (u)ux =0, 

Direct calculation of the left-hand side of the above equation and substitution u t — —ip(u)u x result in 
an expression containing the first partial derivatives of ^\ and r\ and also the factors u x and u x . Since 
the resulting expressions must equal zero for all values of u x , the terms with different powers of u x must 
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be equated to zero separately. Having grouped the terms by powers of u x , one notices that the terms with 
u x cancel each other. Equating the terms without u x and the terms with u x to zero leads to the following 
system of linear PDEs for the unknown components of X: 

rj t + ip{u)r/ x = 0, (3) 
<p(u)^-^-<p(u)^+^(u)^+ V /(u)r, = 0. (4) 

Having solved the above system, one obtains the following infinitesimal operator describing all o.p.g. admit- 
ted by the nonlinear advection equation: 

X = £ (t) d t + (A + <p{u){W + tcp'(u)B) d x + B d u , (5) 

where £w = £^(t, x, u), A — A(u, x — <p(u)t) and B = B(u, x — <p(u)i) are arbitrary smooth functions of the 
indicated arguments. 



2.2 Finite difference discretization and differential approximation 

A discretization of ([j]) is performed on a stencil composed of three (non-collinear) nodes of an unstructured 
grid - (to,Xo), [ti,%i) and (t2,X2). Denoting the values of the grid function at the nodes by Ui and also 
= Xi — Xo, (i = 1, 2), where e is a small parameter, the finite difference discretization can 

be written as: 

T ,. , j \ h%{u\ - wo) - h\{u% - u ) ~.t 2 {u 1 -uq)-t 1 (u 2 -uq) 

L h (to,xo,uo,eT 1 ,ehi,ui,eT2,sh2,u 2 ) = ; ; h (p[u) = 0. (6) 

T\h2 - Tih\ Tih\ — T\h<x 

One can eliminate u\ and u 2 from (^) using Taylor series expansions centered at the node (toj^o)- Omitting 
the index in the resulting expression, one can write the result as 

L +sLi +o(e) = 0, (7) 

where 

L Q = L(x,t,u,u t ,u x ) = u t + tp{u)u x , 

h 2 - ip{u)r 2 2 2 hi - tp(u)Ti 2 2 

M = 7Tr — : r^\ T i u tt + 2T~ihiu tx + h±u xx ) - — — T-^\ T 2 u tt + 2T 2 h 2 u tx + h 2 u xx ). 

l[T\h 2 - T 2 h\) l\T\hi — r-xh\) 

Equation (0) is called a first differential approximation 0] or an equivalent differential equation H , and it 
can be interpreted as a partial differential equation with a small parameter which models the perturbations 
of the original advection process due to the effects of finite difference discretization. Along with the small 
parameter e, equation (Q) has the parameters Ti and hi (i — 1, 2), which are not necessarily small. 



2.3 Approximate Group Analysis of the First Differential Approximation 

Following JIo|, [iljl , consider a small perturbation of the infinitesimal operator X = X$ + eX\ + o(e). This 
perturbed operator corresponds to an approximate o.p.g. (a.o.p.g.) of transformations z' = fo(z, a) + 
efi(z,a) + o(e) by virtue of an "approximate" Lie theorem |10(| . The main problem of approximate group 
analysis is to find all a.o.p.g. of transformations such that the equation ([?]) remains invariant up to o(e) terms 
under these transformations. As in the case of "exact" group analysis, this problem can be solved through 
the approximate infinitesimal criterion of invariance (l^, Theorem 3, p. 98] which, for this case, reduces to 
the solution of the problem 



X^'+eXY 1 ){L + eL x ) = o(e), (8) 

/ Lo+eii=o(e) 

where 

Af ) = + & ] d x + m d u + cfd ut + ci x) d Usc 



i = 0,1, and the terms Q ,Q and Q are obtained from the standard second prolongation formulae 
, Theorem 2.36]. As the above expression for X\ ' shows, in addition to prolonging the action of these 
operators to the second derivatives participating in (Q), we must also prolong it to the parameters Tx,T2, hi 
and hi (the grid step sizes), because the latter are transformed according to the transformations of the 
independent variables. Consider, for instance, T\. Since eri = ti — to, we get: 

e£ (Tl) =e(^ ri) +et[ Tl) + o(s)) = ^ t \h,x 1 ,u l ) -^(t ,x ,u ) =<nA(f +/n^ t) )+"(4 (9) 
Equating the terms of the first order of e and performing the same calculations for t 2 , /ii and /12, one obtains: 

d hl) = + hM^, ^ = t 2 d4 x) + h 2 D x ^\ 

Using expansions up to the second order in e in @ one can obtain the expressions for , £x '£i 
and £^ 2 '' in a similar way. 

Substituting the above relations into (||) and using (0) to eliminate Ut — —(p(u)u x — eL\ + o(e), the 
left-hand side of (||) can be rewritten as 

(X L + e(X L 1 + X 1 L )) = m + <p(u) VOx + u x (<p(u)$ - - ¥>(u)£&> + V 2 (u)^l + <p'(u) m ) 

+eL 1 (u x ^ ) - <p(u)u x ^ - mu + $ + + eXoLi + eXxLo + o(e). (10) 

Equating the leading terms in the above expression to zero, as required by (jq), and noticing that the 

(t) (x) 

resulting equation must hold for all u Xl one obtains that the functions £o j£o an d Vo must satisfy the 
system of equations (§)-(§), resulting from application of infinitesimal criterion of invariance to the original 
equation (EI). 

Continuing decomposition by powers of e in (|l0|), consider the equation obtained by equating the 0(e) 
terms to zero. First notice that one can eliminate u t from the resulting equation using the substitution 
Ut — —ip(u)u x , because taking the term eL\ from (Q) into account contributes only to the terms 0(e 2 ) 
which are not considered here. After u t is eliminated , one is left with a partial differential equation for 
the unknown functions jCi^ an d Vi which depend on independent variables t, x and u. The equation 
also contains r,-,/ij, (i = 1,2) and u x ,Uu,Ut x ,u xx . Since (|^) has already been used to eliminate Ut, the 
resulting equation must hold for all values of the parameters Ti,hi, (i = 1,2), u x ,utt,ut x ,u xx . This allows 
one to decompose the equation further. To this end, notice that among all 0(e) terms in (f_0|), only the term 
eXiLq is independent of Ti, hi, (i — 1, 2). Therefore, the equation can be decomposed to obtain the following 
system: 

AiLq = 0, (11) 
L^u^-^u^-^+^+^u^+XoLx = 0, (12) 

where u t is replaced by — ip(u)u x . It follows from equation ( |TT| ) that in order for the a.o.p.g. Xq + eX\ +o(e) 
to leave equation (j?j) invariant up to 0(e) terms, it is necessary that the components of X\ satisfy the same 
system of equations ([j])- (|]) as the system that an infinitesimal operator of an o.p.g. admitted by advection 
equation (p satisfies. It was also established above that Xq must be admitted by the original equation 
(0). Notice that @ imposes additional conditions on X$. Knowledge of these conditions allows one to 
select those of the o.p.g. (ph of the advection equation which can be made a.o.p.g. of the first differential 
approximation at the price of adding a small correction eX\ to the group in order to account for the effects 
of discretization. Following the terminology of jnj, these groups are said to be inherited by a perturbed 
equation (0). Notice that an infinitesimal operator X^ with all zero components satisfies (|Tl|). Therefore, 
an operator which is inherited by (Q) leaves this equation invariant up to 0(e) terms even without a small 
correction. The purpose of the following is to derive this set of conditions on Xq and to describe all groups 
inherited by (0). 

Substituting the expression for L\ into the left-hand side of (|lj) and multiplying both sides of the 
resulting equation by 2(T\h 2 — T 2 h\), one obtains an equation with zero on the right-hand side and a number 
of terms on the left-hand side. Each of these terms has a factor of the form TiT^h^h^, where k, I, m and n 
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can be 0, 1 or 2, and the terms with the same factors can be grouped together. Since ( |l2| ) must hold for all 
values of ti,T2, hi and h 2 , each group of terms must be equated to zero independently. This procedure gives 
the following system of equations: 

Co" } + 2u tt D t ^ ) + 2u tx D t Q x > - ip{u)u x u tt S^l + u x u tt £^ - UuVou = 0, 
-9f('")(Co* t) + 2tift ACo + 2u tx D t ^ x) ) + ip{u)utt£ i [ x ) + UttCot _ - u tt ip'(u)r) 

+tp 2 (u)u x utt£ ( u + (p(u)u tt r]o u - ip(u)u tt €ot - <fi 2 (u)u tt ^ = 0, 

Q + U tx ACo + WttACo + UxxDtZo^ + UtxDxZo - P(u)u x Ut x £o£ + u xUtx^ u - u txV0u = 0, 

-VJ(u)(C^ te) + UteAf^ + w«A^* } + u xx D t ^ x) + u tx D x ^ ] ) + ^(ujufa,^ + UtxCo? 
-ip(u)u x u tx ^ ~ tp'{u)u tx i] + ip 2 {u)u x u tx S^l + ip(u)u tx r] 0u - <p(u)u t x£ot - ¥> 2 ( u ) u tx£o x = 0, 

Co + 2lf te ACo + 2 «z:£ ACo^ - V'WUs'UxxCol + U x U x X^qJ - U xx T] 0u = 0, 

-^(u)(Co a!x) + 2utxACo + 2u xx D x Q x >) + tp(u)u X x£<te + u xx£,of ~ V{ u ) u xu xx ^ - u xx ip'(u)ri 

+(p 2 («)%«»4t + PM^bIOu - <P(u)Uxx€ot ~ <P 2 ( u ) u xx€ox = 0- 

After the expressions for Co""* , Q and Co an< ^ the total derivatives occuring in the above equations are 
substituted, and Ut is replaced by —ip(u)u x , one finds that in each equation the terms containing the second 
partial derivatives of u cancel each other. Proceeding as before and decomposing each of the above equations 
results in the following set of necessary and sufficient conditions, under which Xq is inherited by (Q): 

Vott = 0, (13) 

Votx = 0, (14) 

Voxx = 0, (15) 

2ip{u) mtu -ip{u)^ t + ^ = 0, (16) 

VOtu - Cote - V(u)V0xu + f(u)^tx = °J ( 17 ) 

IVOxu- tit +<p(u)$i = 0, (18) 

- M»)$L + <p(u)Vo«u = 0, (19) 

-^\u)dt-^l-^)muu + ^(u)^ u + ^(u)^ x J u = 0, (20) 

M^L-^ol + Vouu = 0, (21) 

= 0- ( 22 ) 

Since Co* , Co^ and 770 satisfy (||)-(f|), one can write Co = €0 (t,x,u), Co*^ = A(u, x— ip(u)t)+ip(u)^ (t, x, u) + 
tip' (u)B(u, x — ip(u)t) and 770 = B(u, x — ip(u)t), where C , A and B are arbitrary smooth functions of their 
arguments. Introducing new independent variables x\ — it, 22 = x — ip(u)t,xs = t and denoting the partial 
derivatives of the above arbitrary functions by subindices (eg. A\ 2 — A XlX2 etc), system ( [l3|) - (|22| ) can be 
rewritten as: 

B22 = 0, (23) 
2tpB 12 + 4tp' B 2 - ipA 22 = 0, (24) 
-2<pB 12 -2cp'B 2 + cpA 22 = 0, (25) 
A 22 -2B 12 = 0, (26) 

-2ip'A 2 - 2ipA 12 + 2tipip'A 22 + 2/C03 ~ 2< w 'Co2 + 2<p'' 'B 

+2ip'B 1 - (U(ip') 2 + 3tiptp")B 2 - tttptp'B 12 + tpB u = 0, (27) 
2<pA 12 + tp'A 2 - 2tip<p'A 22 - tp'£$ + 2ipip'^ - tp"B 

-tffBt + (2%') 2 + 3t<p<p")B 2 + At w 'B 12 - tpB u = 0, (28) 
2A 12 - 2tip'A 22 + 2<p'Co2 + W&2 + Uip'B 12 - B n = 0, (29) 
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A n - 2bf/A 12 - t<p"A 2 + t 2 { v ') 2 A 22 + + 2y/$ - 2t(<(/) 2 $ 

+tcp"'B + 2tip"B 1 - M 2 tp'tp"B 2 + tip'B n - 2t 2 (tp') 2 B 12 = 0, (30) 

It follows from (^3) and (p5|) that <p' B 2 = 0. According to this, consider two possibilities: tp' = or £?2 = 0. 

2.3.1 Case ip' = (linear advection) 

Denote y(zt) = c = const. Using (p3h , write _B = B(°^(ari) + i^ 1 ) (a;i)x2, where fl'°)(a;i) and B^(xi) are 
unknown functions. Then (]24|), ([26[), ( p5| ) give A = (B^ (xi))'^ + g'{x\)x 2 + l{x\) with g(zi) and l(x\) 
being unknown functions. Substitution of the above expressions into the left hand side of (|29|) results in 
an expression which is linear in x 2 with its coefficients being functions of x\. This expression can be zero 
only if its coefficients are zero. This results in two equations (B^( Xl ) - 2g(x 1 ))" = and (fi' 1 )^))" = 0, 
and therefore B (0) (i:i) = 2g(xx) + Cx\ + D , B^Jxi) = Exi + F and A = Ex\ + {g{x 1 ))'x 2 + l{x x ) where 
C. D. E and F are arbitrary constants. Finally, ( p0| ) implies that An — 0, and after substituting the above 
expression for A, one obtains (g(xi))'" = and (l(xi))" = 0. Therefore g(xi) = Mx\ + Nx\ + P and 
l(xi) — Kxi + L, with arbitrary constants K, L, M, N and P. 

Collecting the above yields the following expressions for A and B: 

A = Ex\ + (2M + N)x 2 + Kxi + L, 

B = 2Mx 2 1 + (2N + C)x 1 +D + (Ex 1 +F)x 2 . 

Notice that, as follows from ( p3| ) - (|3^), £g can still be an arbitrary function. 

Since the expressions for A and B are linear with respect to eight arbitrary constants, the set of infinites- 
imal operators inherited by (Q) is spanned by the following eight infinitesimal operators 

= ud u , -X"(5) 

X( 3 ) = (x - ct) 2 d x + u(x - ct)d u , X( 7 ) 
X (4 ) = (a; - ct)d u , X {8) 

and an infinite-dimensional set of infinitesimal operators 

X (oo) =$dt+a$>d x , (32) 
where t^p = t^p (t, x, u) is an arbitrary function. 

2.3.2 Case <p> ^ 0, B 2 = 

In this case B = B {0 \xi), where £ (0) (xi) is an unknown function, and (^) immediately gives A 22 = 0, 
and therefore A = A^°'(xi) + x 2 A^ x \xi). Substituting the expressions for A and B into ( 
x 2 ((B^)" — 2(A^Yy (2ip') + f(xi, x 3 ), where j(x\,x 2 ) is an unknown function. After this expression for ^ 
is substituted into (|27|), one obtains — ip' A^ +(p' f 3 + tp" B^ +tp' (B^)' — 0. All terms in this equation, except 
the second one, depend solely on x\. As a consequence, we must have f 3 = g(x%) and / = g{x\)x 3 + l(x\) 
where g{x\) and l{x\) are unknown functions. Using this result, one can write A^ = g + aB^ + (B^)' , 
where a = if" /ip' . 

Collecting together the above results, one obtains 

A^ + (g + aB^ + (B^y)x 2 , 
G{xi)x 2 + g(xi)x 3 + l(xi), 

where G{x\) = —2g — 2{aB^)' — (B^)". Proceeding with the above, substitute the resulting expressions for 

A, B and ffl into ( p8| ) and (|3C|). The first equation will be satisfied identically, while the second will result 
in an equation which has zero on its right-hand side and a linear expression with respect to the independent 



ud x , 
d x , 

u(x — ct)d x + (x — ct) 2 d Ul 
[x — ct)d x + 2ud u 



(31) 



29) gives Co 



i 
so 



.4 
B 

(*) 
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variables xi and X3 with its coefficients being functions of x\ only, on its right-hand side. Equating these 
coefficients to zero results in the following conditions on the unknown functions A^ 3 ' , A^\ B^°> , g and I: 

-<p"A^ +<p"g + 2v'g' + V "'B^ + 2<p"(B^)' = 0, (33) 
(A^)" + <p"G + 2<p'G' = 0, (34) 
(A<°>)" + cp"l + 2<p'l' = 0. (35) 

Substitution of the expressions for A^ and G into ( p3| ) and ([34]) results in a system of two differential 
equations with two unknown functions, B^ and g. Integrating this pair of equations yields 

aB^ 0) +2g = &, (36) 

biB® = C 2 (37) 

where C\ and Ci are arbitrary constants and b = a 2 /2 — a' . Equation (|35| ) can be integrated to find I in 
terms of A^ ' : 

where C3 is an arbitrary constant. In order to formulate the final result, consider two cases: 6 = and 6^0. 
Case 6 = 0. 

In this case ip = K(L — u)^ 1 + M, where L and M are arbitrary constants and if is a non-zero constant, 
since in the present case d + 0. As follows from ©, B^ = B^(xi) can be an arbitrary function and 
C2 = 0. Using ( |36|) and (pq); one finds that the linear space of the infinitesimal operators inherited by (0) 
in this case is spanned by the following infinitesimal operators: 

X (1) = td t +xd x , (39) 

X(2) = dt + <p(<p')v d x , (40) 

*(3) = i8t + ({x-<pt)A<V+<pf)d Xi (41) 

+(-aS^ + (B (0) )' - ^-((aB(°))' + (B^)")^ - <pt) ~ -aB^t V + tip'B^) 8 X 
2 2<p' 2 

+B® d u (42) 

where I is related to A^ by j3^ ) with C3 = 0. Unlike and -X"(2)> the operators X/ s \ and X( 4 ) each 
represent an infinite family of operators, parametrized by two arbitrary functions B^ — _B^ '(u) and A^ = 

Case b £ 0. 

In this case, ([53]) and (Q) immediately give = C* 2 &-2 and 3 = |Ci - ±C 2 a&-3. The functions I and 
as in the previous case, are related through ( pa) . Using these expressions, one can see that, as in the 
cases above, all infinitesimal operators inherited by (hi) form a linear space. For the given case, the basis in 
this space can be chosen as follows: 

X {1) = td t +xd x , (43) 
X {2) = {-^{Aa'b 2 -2ab + 3(b') 2 - 2bb")(x - ipt) - -ab~h) d t 

+(t<p'b-i - hipab-i + [\b-i {a" - 2aa' + 2_) - -^-l(4a'6 2 - 2ab + 3(b') 2 - 2bb")){x - tpb)) 8 X 

2 2 2 8<£> 

(44) 

A (3 ) = {tp')^ d t +<p(<p')i d x , (45) 
X (4) = ld t + {{x-yt)A^ + <pl) d x , (46) 

where I and A^ ^ are related through ( pS| ) with C3 = 0. Unlike the case 6 = 0, only X(4) is an infinite family 
parametrized by an arbitrary function A^ . 
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3 Results and Conclusion 



One parameter groups of transformations approximately admitted by a first differential approximation cor- 
responding to a finite difference approximation of a nonlinear advection equation have been studied. 

It was demonstrated that the necessary condition for such a group to be approximately admitted by a 
first differential approximation is that this group is admitted (exactly) by the original advection equation. 
Sufficient conditions (p"3|)-(|2^), imposing additional restrictions upon the groups of the original equation, 
were also derived. 

The description of a family of transformation groups admitted by a nonlinear advection equation (^) 
contains three arbitrary functions. The requirement that the group must be inherited by a first differential 
approximation (Q) narrows this arbitrariness in the following way. 

For the most general case ip' ^ 0, b ^ 0, the family of transformation groups inherited by (m) is linear, and 
its basis may be chosen to consist of three groups (|43|)-(f45|) and an infinite family of groups (^6|) parametrized 
by an arbitrary function of u only. 

For the case b = 0, the family of inherited groups expands and can be described as a linear set spanned by 
two groups (|39|)-(^0|) and two infinite families of groups (Ji~t|)-(fi2|) each parametrized by a different arbitrary 
function of u only. 

Finally, in the linear advection case when ip' = 0, the set of inherited transformation groups has as its 
basis eight groups (|l|) and one infinite family of groups (|3^) parametrized by an arbitrary function of t, x 
and u. 
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